Neutralino decay rates with explicit R-parity violation 
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In the minimal supersymmetric standard model [Q, 
a discrete symmetry called R-parity is invoked to for- 
bid gauge invariant lepton and baryon number violat- 
ing operators. The R-parity of a particle is given by 
Rp — (_i)i'+2S-i-3B^ where L and B are the lepton and 
baryon numbers, and S is the spin. Standard model 
fermions, Higgs bosons, and gauge bosons have Rp = +1, 
while their superpartners have Rp = — 1. This symme- 
try guarantees the stability of the lightest superpartner 
(LSP). 

There is no deep theoretical motivation for imposing 
R-parity, and it is an interesting exercise to explore the 
phenomenology of R-parity violation We introduce 
explicit R-parity violation by adding 



VK^^ = XLLE" 



X'LQD" + X"U''D''D'^ 



(1) 



to the superpotential. The first two terms violate lepton 
number, and the third violates baryon number. The LSP 
can now decay into standard model particles. 

In models where supersymmetry is broken by super- 
gravity, the LSP is usually the lightest of the neutralinos, 
which are superpositions of the superpartners of the neu- 
tral electroweak gauge bosons and the superpartners of 
the neutral Higgs bosons. R-parity violation allows the 
neutralino, which is a Majorana fermion, to decay into 
three standard model fermions (see figure 1). 

Neutralino decays are relevant to accelerator searches 
for R-parity violation, especially resonant squark produc- 
tion . Astrophysical neutralino decays are also of in- 
terest, and can put strong constraints on the R-parity 
violating couplings [|| . 

It is interesting to note that neutralino decays may be 
relevant to the reported high anomaly at HERA 
The anomaly is an excess of events with a positron in the 
final state at high Q^. One interpretation is resonant 
production of a squark u in e+d u due to an LQD^ 
term in the superpotential. The squark may decay back 
to a positron via the R-parity violating operator or may 
decay into 'x^u and the neutralino then decay into a 
positron by an R-parity violating interaction. This sce- 
nario awaits confirmation, such as from related charged 
current events 

The calculations of neutralino decay rates into 
fermions are subtle because they involve both Majorana 
fermions and fermion-number violating operators. To our 
knowledge, only one calculation is available in the litera- 
ture [Q, and it neglects sfermion mixing. We improve this 
calculation by including the complete sfermion mixings 
and we find some small but significant differences. 

The differential decay rate of the neutralino is given 
by a standard three-body phase space factor multiplied 
by a squared amplitude which is averaged over the initial 
neutralino spin and summed over the final fermion spins. 
We denote this spin sum by a primed summation symbol. 



dxidx2 2567r3 
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Here Xi^2 = 2£'i^2/tox ^^"^ -^i-2 ^''"S two final state 
fermion energies. 

In the following we present the calculation of the spin 
averaged squared amplitudes for neutralino decay. We 
first consider the U'^D'^D'^ term in the superpotential, 
followed by LQD" and finally LLE". 

To calculate the U'^D'^D'^ decay, we first obtain the 
Feynman rule for the R-parity violating vertex. Writing 
all indices explicitly the superpotential reads 



(3) 



where C/f„, D^.^ 



and D'f.^ are the superfields of the right- 
handed quarks (and squarks) respectively, the super- 
script c denotes charge conjugation, i, j, and k are gener- 
ation indices and a, (3, and 7 are SU(3)c triplet indices. 
It follows from the antisymmetry of e"^''' that A"^^, is 
antisymmetric in the last two indices, thus there are 9 
couplings for three generations. The superpotential can 
be written as 



(4) 



with k > j. The Lagrangian that follows from the above 
superpotential contains 



C-UDD 3 



h.c. 



(5) 



where ry^ is the right-handed component of the down- 
quark with generation index i and color index a, and 
dnia is the squark associated to it. To pass to four- 
component notation we use the relation ^1^2 ~ '0i^'i?'02, 
with Pr = i(l -h 75). We obtain 



^UDD 3 



~2e'^^'<X'r, 



URiad'^^pPjidk'^ 



+dRjfju'^^PRdkj + dRkyU'^^PRdjfj + h.c 



FIG. 1. Typical dia- 
gram for neutralino de- 
cay into standard model 
fermions via sfermion ex- 
change. The R-parity vi- 
olating vertex is circled. 
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The sfermion mass eigenstates /«, with k = 1, ... ,6, 
are related to left- and right-handed sfermions through 
the 6x3 mixing matrices tL, ■ and r{ , -, where i = 1, 2, 3 
is a generation index. Without mixing, k = i is a left- 
handed sfermion with generation i and k = i -I- 3 is a 
right handed sfermion with generation i. In general, the 
sfermion mixings are given by 



f 



(7) 



2 



with inverse 
fm 



where 



p/* f 



fu 



p/* f 



(8) 



Note that there is only a left-handed sncutrino, so there 
is only a left handed mixing matrix and k = 1,2, 3. In 
the simple case of no flavor mixing but with left-right 
mixing, the only non-zero elements are 



^ Lii 



cos 6f, 

rL = -ri(,+3)i=sin^/' 



(9) 



where ^/ is the left-right mixing angle for sfermion /. 
Inserting the mixing matrices into the previous La- 

grangian, we have 



+^t.jd.pu<i^PRdk^ + Tt.kd.^ul^PRdip\ + h.c. (10) 

The R-parity violating vertices are read directly from this 
Lagrangian. 

We also need the neutralino-sfermion-fermion vertex. 
It comes from the interaction term 

C^fj- 3 x{9^f,.PL + 9^f,.PR)fiafL + h.c., (11) 

where i is a fermion generation index, k specifies the 
sfermion mass eigenstate, and a is a color index. The 
couplings are 



7^ 



-RA 



R^yxfi 

where A can be L or i? and 



r 



LA 



9xh = - [(9/ - T^) g'Nxi + TzgN^2] 



LR 



RL 
^XUi 



\f2mw sin /3 ' 

\f2m\fi/ cos /3 ' 



LR 
9xdi 



„RL 
9xdi 



(12) 

(13) 
(14) 

(15) 

(16) 



Here Ts is the third component of the weak isospin, Nij 
is the 4x4 neutralino mixing matrix in the convention in 
which all neutralino masses are positive, u is an up-type 
quark or neutrino (T3 = +1/2), and d is a down- type 
quark or charged lepton (T3 = —1/2). The charges are 
Qu = 2/3, Qd = —1/3, = 0, and qe = —1. Notice that 
the only surviving neutrino couplings are . 

We can now obtain the full expression for the decay 
amplitude. In order to get correct signs for the interfer- 
ence terms we use Wick's theorem. The effective operator 
for neutralino decay is 



(17) 



X{G^,Pl + G^^PR)Uic.d'^.pPRdk.y + 

x{G^^Pl + G2fPR)dj0U^PRdk^ + 
x{G^^Pl + Gl^PR)dk^T;^PRdj0 



(18) 



A and B are either i? or L and A^^ = {p^ — m'j + 

im^T^^)~^ is the sfermion propagator. In the case of no 
sfermion mixing, ^/ = 0, we find 



(19) 



In all cases G* corresponds to a fermion and G corre- 
sponds to an antifermion in the decay amplitudes. The 
rates for the charge conjugated decays are found by com- 
plex conjugating G and G* everywhere in our expressions 
for the spin summed squared amplitudes and thus are 
identical. 

Suppressing color wavefunctions, we obtain the follow- 
ing amplitude for the decay % — > Uidjdk- 



M{x Urdjdk) = {uidjdk\T\x) (20) 

VxiG^^PL + G^^PR)VuJid,PRVd^ - 

^x{G^^Pl + G^^PR)vd,Uu,PRVi, + 
Vx{GZPl + G^^PR)vdJiu,PRVd^ , 

where Uf and Vf are the usual particle and antiparticle 
Dirac spinors associated with quark /. Notice the mi- 
nus sign in the second term in brackets, which comes 
from Dirac statistics. It is comforting that the same rel- 
ative sign is obtained from the familiar Feynman rule 
of the sign of the fermion permutations; the first term 
has the fermions in the order {x,Ui,dj,dk), the second 
in the order (xjdj,Ui,dk), and the third in the order 
(XjdkiUiidj). The second and the third are odd and 
even permutations of the first. This confirms the relative 
minus sign for the second term. 

Squaring the amplitude, averaging over the two spin 
states of the neutralino, and summing over the final an- 
tiquark spin states, we obtain 



Y,'\Mx ^ Uidjdk)\^ = 8c/|A^;.fc|2 X Re { 



(21) 



drdk[i\G^ 



RL\2 



\G 



RR\2' 



X-u + 2G«^G««*m;,m„ 



(IG^^I^ + |G.7r)x • dj + 2G^-G^f*m^m,, 



+u ■ dk 

+u- dj^yi'^^k I -^I'^dk I )X- fJ-k-r ^'^d^ 
-G^''Gff*g{u, dk, dj, x) - G^''Gff*idk ■ dj)m^mu 
-G^RG^^*{dk ■ u)m^md, - G^^G^^*(x • 4)m„md, 
_GRRGZ*g{u, dj, dk, x) - G^'^G^^^idk ■ dj)m^mu 
-GTgI'^^u . dj)m^md, - Gf ^G«f *(X • rfi)m„md, 
-Gl^Gl^*g{dj,u,dk,x) - G^^Gl^*{dk ■ u)m^ma, 
-G^^G^^*{dj ■ u)m^md, - G^^Gf/'*(x • t.)m,,.m,, }, 
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with the color factor c/ = 6, and g{a, b, c, d) = [a ■ b){c ■ 
d) — {a ■ c){b ■ d) + {a ■ d){b ■ c). Particle four- moment a 
have been denoted by the particle letter and unambigu- 
ous indices have been suppressed. 

The calculation of the LQD'^ decays is very similar. 
The superpotential is 



WlQD = e-^P [KjkLiaQjpaDl^] 



(22) 



where a and p are SU(2)i, indices and a is an SU(3)c 
index. Suppressing color wavefunctions, we obtain the 
Lagrangian 



^LQD 9 Kjk eudkPLUj + ULjdkPLSi + d-RketPLUj 



-VLidkPhdj - dLjdkPhVi - d}ij,iyfPLdj 



+ h.c. 



(23) 



We have used the identities ^1772 = ■02-PlV'i and ^1^2 = 
'iJjIPliI'2 to pass to four-component notation. Here there 
are 27 couplings for three generations as the coupling 
matrix is unconstrained by symmetry arguments. 

After introducing sfermion mass eigenstates through 
Eq. (H), multiplying by Eq. (|Tl|), and using Wick's the- 
orem, we obtain the following amplitudes for the decays 
X efujdk and x ^ v^d^dk- 



M{x^e+u,dk)= (24) 

V^iCi^^PL + G^,I^PR)VeJld,PLVu, - 



■^ijk 



VxiGt.PL + Gt;'PB)Vu,Ud,PLVe^ + 
UdAGl''*PL + G"^^*Pn)u^Ue.PLV^, 



M{x^ Vid-jdk)-- 



(25) 



A' 



ijk 



v^Gl!" PhV^JidPLV^ ~ 

MGdfPL + G'ifPR)vaMa,PLVu, + 



nddGZ*PL + G^;:*PR)u^u.^PLVd, 
The matrix elements squared then follow as 

Y.'\m{x - e+u,dk)\ = 2c/|A^,p X Re { (26) 

spins 

e • d [(IG^^p + ICe^Hx • u + 2Gl^Gl^*m^m^\ 

• e [(|G«^*P ^ |Gn')x • d + 2G^^*Grm^m<i] 
+y • d [(IGf'^p + IG^^Hx • e + 2G^^G^^*m^me] 
-G^^Gf «g(zi, e, d, x) ~ G^'^Gf «(e • d)m^m^ 
-G^^Gf^(e • u)m^mi - G^«Gf^(x • e)m„md 
-G^ ^G^^*g(e, d, u, x) - G^^G^^* (e • d)m^m^ 
-G^^Gf'«*(u • d)m^me - G^^Gf'«*(x • d)m„me 
-G««*Gf'^*<7(e, u, d, x) - G«^*G,"*(e • u)m^ma 

-G««*Gf'«*(d • u)m^m, - G«^*G,^«*(x • u)mam,], 



2c/|A^,pxRe{ 



(27) 



V ■ d [(IG^f p -I- IG^fHx • d + 2Gif Gif *m^m,^. 

+d ■ y [(|G«^*p + |G««*nx • d + 2G«^*G^>^m,J 
+H.d|G^^Px-z^ 

-Gi^^G««g(d, V, d, x) - Gif G^«(^. • d)m^m^^ 
-G^df^ti^ ■ d)^x^d, " G«^(x • ^^)m,^mrf, 
-G^^fGt'^*g{v,d,d,x) - G^^G^^*(j. • d)mxm<j, 
-Gjf * G^^*g(j., d, d, x) - Gjf *G^^* (t. • 3)m^m,, } , 

where now the color factor is given by c/ = 3 in both. 
Note that in Eqs. ( p6| ) and (^^ for the LQD'^ processes 
the factor in front is 2, not 8 as in Eq. ( |2l| ) for the 
U'^D'^D'^ processes. This is because the U'^D'^D'^ pro- 
cess is identical for A^^-j, and , whereas each element 
of the LQD'^ matrix A.^^ j, gives a unique channel. 

The calculation of the LLE'^ decay is very similar to 
the LQD'^ case. From the superpotential 

Wlle = e"" [KjkUaLjpEl] , (28) 

we obtain the Lagrangian 

t^LLE 3 '^Kjk [ehiekPLVj + i>LjekPLei + e*j^k(^fPLVj 
-VLiekPLGj - eLjekPhVi - eRi^e^jPLVi] + h.c, (29) 

where j > « by the same antisymmetry argument as in 
the U'^D'^D'^ case. Again there are 9 couplings for three 
generations. 

Passing to sfermion mass eigenstates and using Wick's 
theorem, we obtain the following amplitude for the decay 



Mix- 



4''jC'k ) 



(30) 



VxiGef-PL + G^^^Pr)v,^u,,Plv,^ 

VxG^^^PLVu.Ue^PLVe, + 



ue, [G^^Pl + GI^*Pr)u^u,^Plv,^ 
The matrix element squared is 



E' 

spins 



Mix 



8|A, 



ijk\ 



X Re ■ 



(31) 



e • V [i\G':^\' + IG^ r)x • e + 2G^,-G^«*m^me. 
+e-e|G^^|2x-i. 

+e . . [(|G«^*p + |Gf,^*nx • e + 2Gf,^*Gf>, 
-G,"G^^*g(e, e, x) - Ge^,«G^^*(e • v)m^m,^ 
-G^.^Gl^'gie, v, e, x) - G,^^^Gf^^(x • v)vfi,,m,^ 
-Gf'^^Gf^^(e • j.)m^me, - G^^^G^^(e • v)m^m,. 



'Gt'^Gl^giv.e, e, x) - G^^G«^(i/ • e)m^m,. 



4 



Again, there is a factor of 8 instead of 2 due to the fact 
that Xijk and Xjik allow the same process. Note that if 
i > ji the amplitude gains an overall minus sign, but the 
squared amplitude is identical. 

In the limit of no sfermion mixing, we can compare our 
results with those of Butterworth, Dreiner, and Morawitz 
We differ in several respects. In all three types of 
decays we keep complex conjugations in the couplings, 
which are crucial when the neutralino mass eigenvalue 
is negative, and we do not have the global phase space 
factor of 2(1 - m?/-B?)-i/2 which appears in their rates. 
Our LLE'^ and U'^D'^V decay rates are a factor of four 
larger than theirs, due to the fact that each channel is du- 
plicated in the coupling matrix which doubles the decay 
amplitude. Finally, we find that in the U'^D'^D'^ decay 
their couplings a and h are exchanged, in other words 
gauginos and higgsinos are exchanged, and that in some 
instances, particles and antiparticles have been confused 
(because the a's are the same for particles and antiparti- 
cles but the &'s differ). 

Together with H. Dreiner we have agreed that the 
formulas in this paper are the correct ones. Erratum to 
will be published elsewhere 

In conclusion, we have computed the neutralino de- 
cay rate in R-parity violating extensions to the minimal 
supersymmetric standard model. For the first time we 
include complete mixing among sfermions. Our results 
supersede previously published calculations 

We thank H. Dreiner for taking the time to compare 
with our results and for comments on the manuscript. 
This research was supported by grants from NASA and 
DOE. 
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